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ABSTRACT 

We analyzed all archival data of the low-mass X-ray binary system 4U 1636- 
53 with the Rossi X-ray Timing Explorer (1490 observations). We found a total of 
336 type-I X-ray bursts from this source. From fits to the time-resolved spectra, we 
classified 69 of these bursts as photospheric radius-expansion (PRE) bursts. PRE 
bursts show a characteristic time profile in which the fitted blackbody radius increases 
rapidly at the beginning of the burst, and then drops abruptly close to the peak of the 
burst. The lowest value of the radius after the expansion phase defines the so-called 
touchdown point. We found that in 17 of the PRE bursts, after the touchdown point, 
the blackbody radius increases again quickly after about 1 second, and from then on 
the radius decreases slightly or it remains more or less constant. In the other 52 PRE 
bursts, after touchdown, the radius of the blackbody stays more or less constant for 
~ 2 — 8 seconds, and after that it increases slowly. Interestingly, those PRE bursts 
in which the blackbody radius remains more or less constant for >2 seconds show 
coherent oscillations in the tail of the burst, whereas those PRE in which the blackbody 
radius changes rapidly after touchdown show no coherent oscillations in the tail of the 
burst. From a Kolmogorov-Smirnov test we find that the difference between the two 
groups of PRE bursts is significant at a 5-cr level. This is the first time that the presence 
of burst oscillations in the tail of X-ray bursts is associated with a systematic behaviour 
of the spectral parameters in that phase of the bursts. This result is consistent with 
predictions of models that associate the oscillations in the tail of X-ray bursts with 
the propagation of a cooling wake in the material on the neutron-star surface during 
the decay of the bursts. 

Key words: stars: neutron — X-rays: binaries — X-rays: bursts — stars: individual: 
4U 1636-53 



1 INTRODUCTION 



Thermonuclear, type-I, X-ray bursts (e.g., Lewin et al.lll99j : 
IStrohmaver fc Bildstenll2003l : iGallowav et al.ll2008l ) are due 
to unstable burning of H and He on the surface of accreting 
neutron stars in low-mass X-ray binaries (LMXBs). Some X- 
ray bursts are strong enough to lift up the outer layers of the 
star. During these so- called photospheric radius expansion 
(PRE ) bursts (e.g., iBasinska et all 1 1984 iKuulkers et al.l 
I2OO2D . the radiation flux that emerges from the stellar sur- 
face is limited by the Eddington flux. 
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One of the best studied sources of X-ray bursts is the 
LMXB 4U 1636-53. For instance, from observations with 
the Rossi X-ray Ti ming Explorer (RXTE) up to May 2010, 
I Zhang et all (|201ll ') detected 298 X-ray bursts. Most of these 
bursts have standard, single-peaked, fast rising and expo- 
nentiall y decaying light c urves; 52 of these bursts are PRE 
bursts (zhang et al.li201lf ). 



Some of the bursts in 4U 1636-53 show millisec- 
ond oscillations at 581 Hz, the so-called burst oscillations 
IStrohmaver et al] 1 19981 ). This oscil lations likely reflect the 
spin frequency of the neu tron star jStrohmaver et allll997l : 
IChakrabartv et all |2003| ). Similar burst oscillations have 
been detected in several other low-luminosi ty accreting 
neutron-star systems (for a review see, e.g. iMuno et al.l 
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Figure 1. Colour-colour diagram of all RXTE observations of 4U 1636-53. The grey points represent the colours of the source from all 
available RXTE observations. Each point in this diagram corresponds to 256 s of data. The colours of 4U 1636-53 are normalized to the 
colours of Crab. The blue filled circles represent the colours of the persistent emission of the source at the onset of a PRE X-ray burst 
with tail oscillations. The red open squares indicate the same for PRE bursts without tail oscillations. The position of the source on the 
diagram is parametrized by the length of the black solid curve Sa. 



I2OOII ). e.g. 4U 1728-34, 4U 1608-52, KS 1731-260 and 
Aquila X-1. Burst oscillations do not occur in every burst 
from these LMXBs; but when burst oscillations are present, 
they occur sometimes during the rise, sometimes in the 
decay, and sometimes both in the rise and the decay of 
the burst. In KS 1731-260, oscillations are only found at 
high mass accretion rate, both in the rise and the decay 
of the burst, and all but one o f the bursts with oscilla- 
tions also show radius expansion (jMuno et al.ll200ll l. In 4U 
1728-34, burst oscillations (both in the rise and the de- 
cay) are also only detected when the mass accretion rate 
is high, whereas most PRE bursts occur when the accretion 
rate i s low, and these PRE bursts show no burst oscilla- 
tions ijvan Straaten et al.|[200ll : lFrancdl200lD . In 4U 1636- 
53, the situation is more complex than in 4U 1731-260 and 
4U 1728-34. Burst oscillations in 4U 1636-53 are observed 
both in PRE and non-PRE bursts, and are detected both at 
low and high mass accretion rate jZhang et al.ll201ll ). From 
these results, it appears that in 4U 1636-53 burst oscillations 
are neither correlated with mass accretion rate nor with the 
PRE phenomenon. 

Burst oscillations have been explained as arising 
from rotation of a brightness asymmetry on the neutron- 
star surface at the sp in frequency of the neutron star 
(|Strohmaver et al]| 19971 ). Asymmetries in the emission pat- 
tern of the neutron-star surface in the rising phase of ther- 
monuclear X-ray bursts can be due to initially localized nu- 
clear burning at the place where the burst first ignites; the 
flame front subsequently spreads to the entire neutron-star 



surface, and the asymmetry, and hence the oscillatio n s, dis - 
appears (|Strohmaver et al.lll996l '). IStrohmaver et al.l (| 19971 ) 
found that the amplitude of the burst oscillations in 4U 
1728-34 decreases monotonically as the burst flux increases 
during the rising phase of the burst. This result is consistent 
with the spreading hotspot model, since as the spot grows 
in size, the amplitude of the oscillation should decrease. 

As we already mentioned, oscillations are detected not 
only in the rising, but also at the peak and the decaying 
phase (the so-called tail) of X-ray burst; in fact, burst os- 
cillations are most commonly de tected in the tail of the 
bursts (hereafter tail oscillations; van Straaten et al]|200ll : 
iMuno et ai]|2002l : iGallowav et al.ll2008l 'l. Most burst oscilla- 
tions exhi bit a frequency drift of ^ 1 — 5 Hz in the tail of 
the b urst (|Wiinandsll200ll : lGallowav et al.|[200ll : lMuno et all 
l2002t ). In general the frequency of the oscillations increases 
towards an asymp totic va l ue in the t ail of the bu r st, al - 
though Strohmave r et al.1 l| 19981 ') and IStrohmave3 (|l999l ) 
found that in 4U 1636-53 the frequency of the oscillations 
decreases in some bursts. The spreading hotspot model can 
neither explain the tail oscil l ations, nor this freque ncy drift 
(ICumming fc BildstenI I2OO0I : ICumming et al.' 2002) . 

Regarding the tail oscillations, Payne fc Mclatoa, (|2006l ) 
proposed that during the decaying phase of the burst, the 
burning front is stalled by the presence of a magnetic field; 
the combination of partial surface burning and magnetic 
fields could lead to anisotropic emission during the tail 
of X-ray bursts. Alternatively, a cooling wake in the tail 
of the burst due to hydrodynamic instabilities can also 
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produce asymmetric emission (|Spitkovskv et alj 12002*). Fi- 
nally, instability modes (eg., presure, gravity, buoyancy, 
etc.) excited in the neutron star burning layer can also 
produce burst oscillations; this scenario can also account 
for the observed fr equency drift of the oscillations in the 
tail of some burst sJCumming fc Bildstenll200(]| : lHevlll2004l : 
IPiro fc Bildstcn 20051 ^ 

In this paper we compare simultaneous power density 
spectra and time-resolved energy spectra of 336 X-ray bursts 
in 4U 1636-53. We find that bursts with oscillations in the 
tail of the burst always show an extended period of a more 
or less constant blackbody radius during the burst decay. 
We describe the observations and data analysis in Sj2l and 
we present our results in ^ Finally, in Sj4] we discuss our 
findings and compare them with current models for burst 
oscillations. 



2 OBSERVATION AND DATA ANALYSIS 

We analyzed all available data (1490 observations) of 4U 
1636-53 from the Proportional Counter Array (PCA) on 
board RXTE. The PCA consists of an array of five coUi- 
mated proportional counter units (PCUs) operating in the 
2—60 keV range. We produced 0.5-s light curves from the 
Standard-1 data (0.125-s time-resolution with no energy res- 
olution) and searched for X-ray bursts in these light curves 
following the procedure described in Zhang ct al. ( 201l] ). 

We used the Standard-2 data (16-s time-resolution and 
129 channels covering the full 2-60 keV PCA band) 
to ca lculate X-ray colours of the source (see IZhang et al.l 
I2OIII . for details). Hard and soft colours are defined as the 
9.7 - 16.0/6.0 - 9.7 keV and 3.5 - 6.0/2.0 - 3.5 keV count 
rate ratios, respectively. We show the colour-colour diagram 
(CCD) of all observations of 4U 1636-53 in Figure [T] We 
parametrized the position of the source on t he diagram by 
the l ength of the solid curve Sa (see, e.g. iMendez et al.l 
119991 ). fixing the values of Sa = 1 and Sa = 2 at the top- 
right and the bottom- left vertex of the CCD, respectively. 

In order to study the bursts in detail, we used the high- 
time resolution modes that were available for each observa- 
tion to produce time-resolved spectra of each burst. About 
8% of the observations have a mode with 500-/is time res- 
olution. The rest of the observations have a mode with at 
least 125-/is time resolution. For every burst we produced a 
spectrum every 0.25 s during the rising part and the peak, 
and 0.5 s, 1 s and 2 s during the decay, covering the whole 
duration of the burst. The variable exposure time was tai- 
lored to collect more or less the same number of counts per 
spectrum. We generated the instrument response matrix for 
each spectrum with the standard FTOOLS routine pcarsp, 
and we corrected each spectrum for dead time using the 
methods supplied by the RXTE team. Because of the short 
exposure of each spectrum, in this case the statistical errors 
dominate, and therefore we did not add any systematic error 
to the spectra. For each burst we extracted the spectrum of 
the persistent emission just before (or after) the burst to use 
as background in our fits; this approach, used to obtain the 
net emission of a burst, is a well established p rocedure in X- 
ray burst analysis fe.g. iKuulk ers et al.ll200j ). We note that 
this procedure fails if the blackbody component during the 
burst comes from the same source that produces the black- 



body component seen in the persistent emission, since the 
difference between two blackb ody spectra is not a blackbody 
(|van Paradiis fc Lewin|[l986f ). This effect is significant only 
when the net burst emission is small, and therefore problems 
may arise only at the start and the tail of the burst, when 
the burst emission is comparable to the pe rsist ent emission 
(see t he discussion in lKuulkers et al]l2002l '). In lZhang et al.l 
we already established that this issue does not affect 
the spectral results in 4U 1636-53. 

We fitted the spectra using XSPEC version 12.7.0 
(|ArnaudllT99^ ). restricting the spectral fits to the energy 
range 3.0 — 20.0 keV. We fitted the time-resolved net 
burst spectra with a single-temperature blackbody model 
{bbodyrad in XSPEC), as generally burst spectra are well fit- 
ted by a blackbody (Gallowav ct al. 200^. We also included 
the effect of interstellar absorption along the line of sight 
using the XSPEC model wabs. During the fitting we kept 
the hydrogen colum n density, A^h, fixed at 0.36 x lO^^cm^^ 
JPandel et al.|[200^ ). and to calculate the radius of the emit- 
ting blac kbody area in km, -Rb b, we assumed a distance of 
5.95 kpc (|Gallowav et al.ll2008l ). 

For each burst we computed Fourier power density spec- 
tra (PDS) from 2-s data segments for the duration of the 
burst using the data over the full PCA band pass, setting the 
start time of each segment to 0.125 s after the start time of 
the previous segment. Because of this, the individual power 
spectra are not independent. The Nyquist frequency of these 
power spectra was always 1024 Hz. We used these PDS to 
produce ti me- frequency p lots (also known as dynamic power 
spectra; see lBerger et al.l l996) of each burst. For each burst, 
we searched for coherent oscillations between 577 Hz and 582 
Hz for the duration of the burst. We considered that a sig- 
nal was significant if it had a probability of < 10~^ that it 
was produced by noise, accounting for the number of inde- 
pendent trials. Since the individual power spectra are not 
independent, to estimate the number of trials we divided 
the duration of the burst by the length of the PDS (2 s). 
Finally, we consider that a burst has tail oscillations if the 
oscillations appear after the peak o f the burst. If we no rmal- 
ize the power spectra according to iLeahv et all (|l983l ). the 
fractional rms amplitude at a given frequency is 



(1) 



where is the measured power, Ij is the total number of 
counts (source plus background), and is the number of 
background counts jBelloni fc Hasingeilll99T} ). To calculate 
the signal power, P^, from the measured power, and account- 
ing for the distribution of powers from Poisson noise in the 
power spectrum, w e used the algorithm described in the ap- 
pendix of Vaugh an et al.l lfl994'). We calculated the error 
ba rs of the rms ampli tude using the same method discussed 
in IWatts et al.l l|2005h . For details on how t o calculate the 
signal power from t h e mea s ured power, see Vaughan et al.l 
(19941), iMuno et all (|2002[ ). IWatts et al.l l|2005l ) and IWattsI 
(2013 ). 



3 RESULTS 

We examined the time-resolved energy spectra of all the 
336 bursts in 4U 1636-53 to identify the PRE and non-PRE 
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Obsid 10088-01-08-01 
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(a) Burst with tail oscillation 




10 15 
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(b) Burst without tail oscillation 



Figure 2. Left panel: A PRE burst with tail oscillation. Right panel: A PRE burst without tail oscillations. The black histogram shows 
the average count rate in 2-s interval every 0.125s. The contour lines show constant power values, increasing from 10 to 40 in steps of 10 
(values are in Leahy units), as a function of frequency (see the left y-axis). Black filled circles connected by a line show the blackbody 
radius as a function of time (see the right y-axis). The red vertical lines define the PTD phase (see text). Note also the power contours 
at ~ 579 — 581Hz at the beginning of the bursts, which are due to oscillations in the rising of the burst. 



bursts. Here we concentrate only on the PRE bursts, and 
in particular in the time interval immediately after the ra- 
dius expansion and contraction phase. In all PRE bursts the 
blackbody radius first increases, it then decreases abruptly 
to a local minimum (the so-called touch-down, TD, point) 
and after that it either increases or decreases slowly. We 
show two examples of PRE bursts in Figure [2] In order to 
compare the power spectrum and the fitted blackbody ra- 
dius within the same time interval, we first fitted the energy 
spectrum every 0.25 s, and we then calculated the average 
blackbody radius in 2-s interval (8 consecutive energy spec- 
tra). We plot the average blackbody radius as a function of 
time with black filled circles connected with a line in Figure 
[5] We also show the burst fight curve at a resolution of 2 s 
and a 0.125-s step with the solid line in that Figure. 

We find that the behaviour of the blackbody radius after 
TD is not the same in all PRE bursts. In 52 out of the 69 
PRE bursts, after the expansion phase the blackbody radius 
first decreases rapidly, it then continues decreasing at a lower 
rate, it reaches a minimum value of ~ 7 — 8 km, and finally 
it increases slowly towards the tail of the burst (see Figure 



2(a) I. In the other 17 PRE bursts, after the expansion phase 



the blackbody radius first decreases rapidly to a minimum 
of ~ 7 — 8 km, then it immediately increases again very 
quickly, and finally it either decreases slightly, or it remains 



the post touchdown (PTD) phase as the contiguous time 
interval after the peak of the burst in which the radius of 
the fitted blackbody is less than 8 km. Our results do not 
change significantly if we choose a value between 8 km and 
8.5 km. If we used a value smaller than 8 km, the length 
of the PTD phase of several bursts would be zero, while if 
we used a value larger than 8.5 km, the length of the PTD 
phase of several bursts would be unbound. We indicate the 
PTD phase by the two red vertical lines in Figure [2] We 
find that most PRE bursts in 4U 1636-53 show a long PTD 
phase, fpDT > 2 — 8s. 

We also examined all the dynamic power spectra of 
these PRE bursts, concentrating only on the oscillations in 
the decaying phase of the burst. We find that 52 out of the 
69 PRE bursts in 4U 1636-53 have tail oscillations (see con- 



tour plot in Figure 2(a) \. The contours of constant power 



as a function of frequency and time shown in Figure 2(a) 
were calculated from 2-s power spectra. We calculated the 
upper limit of the power for the 17 PRE bu r sts in which we 
did not detect tail oscillations l|Grothlll975l : IVaugh an et al.l 
1 1994 ). Except for four bursts, the upper limits are lower than 
the average power of the detected tail oscillations. During 
the other four bursts only one or two of the five PCU de- 
tectors were on, and hence the upper limits are not very 
constraining. From the two examples in Figure[5]it is appar- 



more or less constant (see Figure 2(b) l. We can classify all ent that the burst in Figure 2(a) which shows a long {^^ 6s) 



69 PRE bursts in 4U 1636-53 into one of these two groups. 

While it is apparent that the duration of the phase 
around the minimum radius is not always the same among 
the 69 PRE bursts, we need to find an objective way to 
measure the length of this phase, for instance choosing a 
contiguous time interval within which the radius is below a 
certain value. Clearly this value has to be larger than the 
minimum radius reached in all bursts after the expansion 
phase, and it has to be smaller than the local maximum 
of the radius just after the radius at the touchdown point 



PTD phase, has oscillations at the tail, while the burst in 
Figure 2(b) which has a short (~ Is) PTD phase, has no 



oscillations at the tail. 

We find that bursts that show tail oscillations have in 
general longer PTD phase than those that do not show tail 
oscillations. In order to quantify this, we calculated the dura- 
tion of the PTD phase, hereafter the PTD time, for all PRE 
bursts and divided the PRE bursts into two groups; burst 
with and without tail oscillations. Figure 3(a) shows the dis- 



in bursts like the one in Figure 2(b) We therefore define 



tribution of the time for PRE bursts with (blue thick lines) 
and without (green thin lines) tail oscillations. This plot con- 
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Figure 3. Top panel: Distribution of, respectively, the PTD time and average PTD radius for the PRE bursts with and without tail 
oscillations in 4U 1636-53. Bottom panel: Distribution of the Sa values for, respectively, the raw data and exposure-normalized bursts 
with and without tail oscillations in 4U 1636-53. The bursts with and without oscillations are shown by thick blue lines and thin green 
lines, respectively. 



firms our initial impression: bursts with tail oscillations have 
on average ~ 4 times longer PTD times than burst without 
tail oscillations. We carried out a Kolmogorov-Smirnov (KS) 
test to assess whether the two distributions are consistent 
with being samples of the same parent population. We find 
a chance probability of 3.5 x 10~^. 

We also calculated the average PTD blackbody radius 
(hereafter PTD radius) for PRE bursts with and without 
tail oscillations. The two distributions are plotted in Figure 



3(b) The KS-test probability that both samples come from 
the same parent population is 2.2 x 10~^. 



Figure 3(c) shows the distributions of Sa for the PRE 
bursts with and without tail oscillations. We find that the 
PRE bursts with tail oscillations have Sa values that are 
larger than 1.9, and the distribution peaks at Sa ^ 2.1, 
whereas PRE bursts without tail oscillations distribute uni- 
formly from Sa ~ 1.7 to Sa ~ 2.4, 

To compensate for the fact that RXTE did not sample 
the CCD of 4U 1636-53 evenly, we normalized the bursts 



number per Sa bin in Figure 3(c) by the total exposure time 



with RXTE at each position in the CCD. We show the re- 
sulting distibution in Figure |3(d)| We find that the distri- 
bution of Sa in PRE bursts with tail oscillations still peaks 
at Sa ~ 2.1, whereas the distribution of Sa in PRE bursts 
without tail oscillations peaks at Sa ~ 1.75. The KS prob- 
ability that the two Sa distributioni|3 come from the same 
parent population is 4.4 x 10"''. 

In order to check whether there is a correlation between 
the PTD radius and the amplitude of the tail oscillation, 
for each burst we calculated the PTD radius and the rms 
amplitude of the tail oscillations every second and averaged 
the values for each burst. After that, we rebined the data by 
a factor of 6 and plotted them in Figure |4l From this Figure 
it appears that the fractional rms amplitude decreases as the 



' This is the KS probability from the raw data, i.e., without 
normalizing the number of bursts per Sn interval according to the 
RXTE exposure along the CCD. We get an even lower probability 
if we instead compare the two histograms in Figure [3(c) [ using the 
test. Here we take the most conservative result. 
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4 DISCUSSION 
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Figure 4. Fractional rms amplitude versus PTD radius for all 
the PRE bursts with tail oscillations in 4U 1636-53. 
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Figure 5. A non-PRE burst with tail oscillations in 4U 1636-53. 
The symbols are the same as in Figure [S] 



average PTD radius increases. We fitted the data both with 
a constant and a hnear function, and we carried out an F-test 
to compare both fits. The F-test probability is 6.1 x 10"'^, 
indicating that a hnear fit is marginally better than a fit 
with a constant. We found no significant correlation between 
fractional rms amplitude and PTD time. 

We also detected nine non-PRE bursts with tail oscilla- 
tions in our observations. Similar to the case of PRE bursts, 
after the peak of the burst, the energy spectra of these non- 
PRE bursts show a period in which i?bb remains more or 
less constant during the time in which tail oscillations are 
present (see Figure [5}. However, in this case it is difficult to 
identify the PTD phase because non-PRE bursts do not have 
(by definition) a radius expansion phase, and a subsequent 
TD point. We therefore did not include non-PRE bursts in 
our analysis, although it is quite possible that the connec- 
tion between constant Rhh and tail oscillations applies also 
to this kind of bursts. 



We analyzed all 336 type-I X-ray bursts in the LMXB 4U 
1636-53 observed with RXTE; 69 of them are PRE bursts. 
For the first time, we correlated the behaviour of the spectral 
parameters of the bursts with the presence of oscillations in 
the decaying phase of these PRE bursts. We find that, after 
the radius contraction phase, in some bursts the blackbody 
radius reaches a minimum value followed by a fast increase 
(short PTD phase). We do not detect burst oscillations dur- 
ing the decaying phase of these bursts. In other bursts, the 
blackbody radius reaches the minimum value followed by a 
slow evolution (long PTD phase). We do detect tail oscil- 
lations in these bursts. We compare the PTD time in the 
PRE bursts with and without tail oscillations, and find a 
significant difference (5o") between the two burst classes. 

The mechanism that produces burst oscillations, and 
why these oscillations are n ot present in all type-I X-ray 
bursts, still remains unclear ( Strohmaver et al]|l996l . Il998l : 
Muno et al. 2002; Muno 20di). Unstable nuclear burning 
is likely not happening uniformly across the neutron-star 
surface so, as the neutron star rotates, variations of the 
neutron-star surface brightness and the neutron-star ro- 
tation sh ould pro duce oscillations during an X-ray burst, 
[strohmav er et al. I (l996) suggested that burst oscillations 
are caused by asymmetries due to initially localized nu- 
clear burning (the ignition point of the burst) that later 
spreads over the surface of the neutron star in the ris- 
ing phase of the burst. This scenario, however, cannot ex- 
plain the tail oscillations that persist for as long as 5- 
10 s, unless the asymmetry can be maintained for such a 
long period. iBildster] (|l995h suggested that in some bursts 
only part of the surface is burning; magnetic fields and 
patchy burni ng could lead t o anisotropic emiss ion d uring 
X-ray bursts. iCummitig fc Bildsten (2000). Hevli (|2004l ') and 
IPiro &: Bildsten 1 20051 ') proposed that modes (eg., presure, 
gravity, buoyancy, etc.) excited in the neutron-star burn- 
ing layer can produce burst oscillations in the tail of the 
burst, and that the frequency drift observed in the tail of 
some bursts is due to a changing mode fr equency as the 
burning layer cools. ISpitkovskv et all (|2002l ) suggested that 
the oscillations in the rising phase of the burst are from a 
slow-moving spot on the neutron-star surface. The bright- 
ness anisotropy may ar ise from hydrodynamic instabilities 
in the burning material. ISpitkovskv et all (|2002l ) found that 
the speed of the burning front near the equator is higher than 
that near the poles. Furthermore, these authors also sug- 
gested that tail oscillations could be due to the spread of a 
cooling wake, which is formed by vortices during the cooling 
of the neutron-star atmosphere. In this scenario, the speed 
of the cooling wake would also depend on latitude. ICummina 
(200^) connected the ideas of unstable zonal flows and os- 
cillation modes, and studied the linear stability of shearing 
zonal flows during the tail of type-I X-ray bursts. He found 
that a differential rotation of 2% between the pole and the 
equator, with the equator spinning faster than the poles, 
is unstable to hydrodynamic shear instabilities. Growth of 
shear instabilities may explain the brightness asymmetry 
that produces the oscillations in the tails of X-ray bursts. 

Our re sults matc h som e of the predictions of the model 
bv ISpitkovskv et all diooj -). if we assume that the bursts 
with tail oscillations are due to a cooling wake starting near 
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the poles, while bursts without tail oscillations are due to a 
cooling wake starting near the equator. According to this 
model, the width and speed of the cooling wake should 
decrease by a factor of ~ 4 as the front propagates from 
the equator to the pole. If the cooling wake starts from 
the equator, then the entire equator belt is likely covered 
very rapidly, and the asymmetry during the cooling, which 
is needed for the burst oscillations, disappears. After the at- 
mosphere contracts to the neutron star surface, the emission 
area changes very quickly due to the high speed of the cool- 
ing wake near the equator. These bursts would then have no 
tail oscillations and a short PTD phase. 

If the cooling wake starts at high latitude, the front 
speed is slower than that in the equator (see Fig. 8 in 
ISpitkovskv et al.l 12002). After the atmosphere contracts to 
the neutron-star surface, the emission area changes slowly. 
The low speed of the cooling wake on the neutron star means 
that the asymmetric emission during the tail of the burst 
lasts longer, and the emission area changes slowly. These 
bursts would then have tail oscillations and a long PTD 
phase. It is interesting that bursts with tail oscillations have 
a PTD time that is about 4 times longer than that of bursts 



without tail osci llations (see Figu re 3(a ) I, consistent with 
the prediction of ISpitkovskv et a l. (2002) 

Following the mod el by Spitko vskv et al.l l|2002l ). 
ICooper fc NaravanI l|200'it l found that in a neutron-star X- 
ray binary system, when the source evolves from low to 
high mass accretion rate, the area where physical conditions 
for unstable burning (where the burst should ignite) moves 
from low to high latitudes on the neutron-star surface. We 
find that PRE bursts with tail oscillations always appear in 
the colour-colour diagram at high S'a values (see histogram 
with thick lines in Figure 3(c)[). 5a is considered to be corre- 
lated to mass accret ion rate l|Hasinger fc van der Klidl 19891 : 
iMendez et all 1 19991 ). This suggests that the cooling wake 
of PRE bursts with tail oscillations starts at high latitude, 
when the source is in the high mass accretion rate state. 

If bursts oscillations are due to a brightness asymme- 
try on the neutron-star surface, the larger the area of the 
asymmetry the smaller the amplitude of the modulation. In 
Figure S] we find marginal evidence that the fractional rms 
amplitude decreases as the average PTD radius increases in 
bursts with tail oscillations. 

We note that the observed changes in the blackbody ra- 
dius in the tail of X-ray bursts could also be due to the effect 
of spectral hardening (colour factor), caus ed by the electron 
scatt e ring in the neutron star atmosphere l|Suleimanov et al] 
I2OIII : IZhang et al.ll201ll ). Since the colour factor depends 
mainly on the ratio L/L^dd, and the PRE bursts in 4U 
1636-53 always show similar flux during the PTD phase, 
the colour factor cannot explain the different dura t ion o f 
the PTD phase in the PRE bursts. iKaptein et al] (|2000l ) 
tried to correct the temperature and radius in X-ray bursts 
for this effect, but changes in the blackbody radius do not 
disappear after applying this correction. 

The reason to study the tail oscillations only in PRE 
bursts is that the phase after the TD point can only be iden- 
tified in PRE bursts (since non-PRE bursts do not have a 
touchdown point). Nevertheless, the relation between PTD 
radius and PTD time on one hand, and tail oscillations on 
the other, may actually extend to non-PRE bursts as well. 
Some non-PRE bursts in 4U 1636-53 show tail oscillations. 



and the blackbody radius stays constant as well during the 
oscillating time (see Figure [5]). This suggests that tail oscil- 
lations are always associated with an emitting area that re- 
mains constant for a while, regardless of the nature (PRE or 
non-PRE) of the bursts. We note that there are instances in 
which the blackbody radius stays constant for a sufficiently 
long period in some non-PRE bursts in 4U 1636-53, whereas 
tail oscillations are not detected in these bursts. This sug- 
gests that the blackbody radius staying constant for a while 
is a necessary but not a sufficient condition for the presence 
of tail oscillations in 4U 1636-53. 

Both positive and negative drift of the frequency 
of burst oscillations have been detected i n 4U 1636- 
53 fetrohmaver et al.lll998l : IStrohmaveij|l999l ). IStrohmaveij 
( 199!^ found that in 4U 1636-53 an episode of a negative fre- 
quency drift was correlated with the appearance in the burst 
of an extended tail of emission with a decay timescale much 
longer than in other bursts from this source. If we assume 
that tail oscillations are from vortices in the neutron-star 
atmosphere (S pitkovsky et al. 20o3). the direction in which 
the vortices drift on the surface of the neutron-star may 
affect the oscillation frequency. When the vortices move to- 
ward the pole, the frequency of oscillations decreases and the 
low-speed cooling wake makes this a burst with an extended 
emission tail. When the vortices move toward the equator, 
the frequency of oscillations increases and the high-speed 
cooling wake makes these bursts decay fast. 

Our analysis shows that tail oscillations in type-I X-ray 
bursts in 4U 1636-53 are always associated with an emitting 
area that remains more or less constant for at least ~ 2 — 8 
s. A similar trend is apparent in another LMXB system, 4U 
1728-34 (Zhang et al. in prep.). In hin dsight, this trend in 
4U 17 28-34 is already visible in Fig. 1 of van Straaten et aU 
l|200ll ) and in 4U 1731-260 in Fig. 5 of lMuno et all (|2000l ). 
although it was then not recognized by those authors, prob- 
ably because of the low number of bursts available at the 
time. 
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